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Abstract The big mass difference between the pion(7r) and rho meson(p) possibly originates from 
the spin-dependent nature of the interactions in the two states since these two states are similar 
except for spin. Both ir and p are quark-antiquark systems which can be treated using the two- 
body Dirac equations (TBDE) of constraint dynamics. This relativistic approach for two-body system 
has the advantage over the non-relativistic treatment in the sense that the spin-dependent nature is 
automatically coming out from the formalism. We employed Dirac's relativistic constraint dynamics to 
describe quark-antiquark systems. Within this formalism, the 16-component Dirac equation is reduced 
to the 4-component 2nd-order differential equation and the radial part of this equation is simply a 
Schr6dinger-type equation with various terms calculated from the basic radial potential. 

We used a modified Richardson potential for quark-antiquark systems which satisfies the conditions 
of confinement and asymptotic freedom. We obtained the wave functions for these two mesons which 
are not singular at short distances. We also found that the cancellation between the Darwin and spin- 
spin interaction terms occurs in the tt mass but not in the p mass and this is the main source of the 
big difference in the two meson masses. 

Keywords pion mass • two body dirac equation • quark-antiquark bound state 



1 Introduction 

The 7r and p are similar states consisting of an u-quark and a J-antiquark except for the spins as n is 
in the spin-singlet state while p is in the triplet. Therefore we can easily attribute any differences to 
the spin-dependent nature of the interaction in the two states. However their mass difference is quite 
large in comparison with the mass differences between other meson pairs such as r\c an d J/& (cc) or 
Do and D'q (cu). In this ground the large mass difference might also be due to the small mass of ir, 
which is explained as a result of the chiral symmetry breaking in QCD. Therefore it will be helpful to 
identify the origin of the small mass of ir and to trace the big mass difference between tt and p. 

Both tt and p are quark-antiquark systems which can be described within non-relativistic frame- 
work. However quarks are light particles and interact strongly; they should be treated relativistically. 
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Furthermore, spin-dependence is purely relativistic in nature and cannot be derived non-relativistically. 
Therefore we need a relativistic description rather than a non-relativistic description, for the quark- 
antiquark systems. Todorov has suggested the two-body relativistic equations for spinless two-body 
systems pQ. Based on the Dirac's constraint formalism and a minimal interaction structure the TBDE 
has been tested and confirmed in classical [5] and quantum field theory [3J. Within this formalism, 
two coupled 16-componcnt Dirac equation with 2 fermion basis is reduced to the single 4-component 
2nd-order differential equation and the radial part of this equation are simply Schrodinger-type with 
various terms calculated from a basic radial potential, which can be phenomenologically adopted. We 
have constructed the two-body Dirac equations (TBDE) [2J [5] for two interacting fermions in a form 
that is as easy to handle as the non-relativistic Schrodinger equation. The brief description of this for- 
malism is summarized in Sec. 2. For the interactions between a quark and an antiquark we adopt some 
phenomenological potential, which will be described in Sec. 3. Final numerical results and discussions 
will be presented in Sec. 4. 



2 The Relativistic Schrodinger-like equation 

Van Alstine and Crater utilized the two body Dirac constraint dynamics for two interacting fermions 
with mass-shell constraint to reduce the relative energy and time degrees of freedom [2j 0]. They 
applied the Pauli reduction and scale transformation and succeeded in reducing the 16 component 
Dirac equation into a 4 component relativistic Schrodinger-type equation, 

m = £lHl ^ £2? S = {P 2 + $v,(<Ti,<T2,P±,A(r), S(r))} = b 2 (w)^, (1) 

where w is the invariant mass and b 2 (w) = — with the relativistic energy e w and the relativistic 
reduced mass m w . <!> w plays a role of interaction in the Schrodinger-type equation ([I]) and includes many 
relativistic terms: Darwin(D), spin-spin (SS), tensor(T), spin-orbit(S'O), etc. For a scalar potential S 
and a time-like vector potential A, <P W can be expressed as follows [5] : 

<P W = 2m w S + S 2 + 2e w A - A 2 + $ D + L ■ (en + <7 2 )& so + a x ■ <t 2 $ ss (2) 
+en ■ fa 2 ■ fh ■ (en + 02)^ sot + (3<n • fer 2 • f — <n ■ 02)§t 
+L • (cti - a 2 )^soD + iL ■ (cn x a 2 )<Psox- 

If S and A are central radial potentials, this equation is reduced further into two-coupled radial 
equations. For it which is in a spin-singlet state, the differential equation becomes 



dr 2 



2m w S + S 2 + 2e w A -A 2 +<P D - 3<£ ss j «o = b 2 v . (3) 



p is in a spin-triplet state, which is a superposition of the S- and D-states. These two states are coupled 
by the following differential equations 

f d 2 1 2\/2 

+ 2m w S + S 2 + 2e w A -A 2 +$ D + <P SS J. u s + -^-{3<P T - 6$sot)u d = b 2 u s , (4) 



{-^2 +^2+ 2m ^ S + 52 + 2e ^A -A 2 +<P D - 6$so + ®ss - 2<£t + 2<?sot| u d (5) 

9 fo 

+ — {Z<P T )u s =b 2 u D . 

Here vq — ifj^{r)/r, and us and up are corresponding quantities of p for the S-state and D-state, 
respectively. Solving these equations are not easy but at least manageable numerically. 
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3 Phenomenological Potential 

To complete the differential equation, we need a basic potentials between a quark and an antiquark. 
Commonly used potential for two quarks are the static Cornell potential, which is composed of a 
static Coulomb-type potential and a linear potential [6]. However this potential does not reflect the 
asymptotic freedom which is one of the important characteristics of the q — q interaction. Richardson[7] 
has suggested the following form for heavy quarkonia, 

~ . . 167T 1 . . 

^ = ~ln q 2 \n{\ + q 2 /A 2 ) ^' 

in momentum space. After Fourier transform to coordinate system, it becomes 

8^r _ 8^4r) 
^ 1 27 27r y ' 

The function f(Ar) needs to approach — l/\n{Ar) as r — > and to approach 1 as r — > oo for the 
potential satisfying both crucial characteristics in QCD, the asymptotic freedom and the confinement. 
We modified this potential using 3 parameters, A, K, and B, 

_ 8nA 2 r 16tt 
W ~ ^27 27rln(ATe 2 + B/(Ar) 2 )' ( ' 

These 3 parameters will be determined to fit the experimental values of meson masses. The one feature 
of this method is that there is no cutoff parameter introduced. Everything is treated non-perturbatively 
and it gives good results for heavy quark meson systems as well as light two quark systems. 

Partitioning the confinement term in V(r) of Eq.Q as an invariant S responsible for the scalar 
potential and the Coulomb-like term as an invariant A responsible for the full four-vector potential, we 
can obtain the specific forms for the various terms in Eq.Q. The readers can refer the Eqs.(A17)-(A20) 
in the Ref.|5] for their explicit forms. 



4 Results and Conclusion 

We fit the invariant masses of 32 mesons and determine 3 parameters in the potential. Their values 
are A = 0.4218 GeV, B = 0.05081, K = 4.198. Here the quark masses used to give the best fit are 
m u — 0.0557 GeV, rrid — 0.0553 GeV. Then using these values we solve the differential equations, 
Eq.([3| for 7r and Eqs.Q and ([5| for p, and obtain their wave functions. They are shown in Fig.l. The 
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Fig. 1 The wave functions of n and p. The solid (black in color) line is tv wave function. The dashed (red in 
color) curve is the S-state wave function of p and the dotted (blue in color) curve is the D-state wave function 
of p. 
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Table 1 The contributions of all terms in Eqs.|3j-([5| . The column < 00 > is for n with S-state only and the 
remaining four columns are the diagonal and off-diagonal terms for p with both S- and D-state wave functions. 
The numbers in () are experimental values. All terms except w are in a unit of GeV 2 and w is in a unit of GeV. 







P 


Terms 


< 00 > 


< ss > 


<DD> 


< SD > 


< DS > 


< -cP/dr* + J(J+l)/r 2 > 


0.8508 


0.3085 


0.2812 






< 2m w S > 


0.0103 


0.0026 


0.0006 






< S 2 > 


0.0942 


0.1631 


0.0446 






< 2e w A > 


-0.0598 


-0.2091 


-0.0211 






< -A 2 > 


-0.4279 


-0.1247 


-0.0094 






Sub Total 


0.3832 


0.1680 


0.0146 






< <Pd > 


-3.804 


-0.2790 


-0.0436 






< -3$ss > 

< &SS > 

< -2<P T > 


3.340 


-0.0364 


-0.0095 
-0.0417 






< -6$so > 

< 2$sot > 
Sub Total 


-0.4643 


-0.3154 


0.1133 
0.0215 
0.0346 






(2V2/3) < 3<P T > 








-0.3088 


-0.3088 


(2V2/3) < -6&SOT > 
Sub Total 








0.6177 
0.3090 


0.3088 


Total Sum 


0.0033 


-0.1750 


0.3158 


0.3090 


-0.3088 


b 2 


0.0033 


0.1411 


w(GeV) 


0.159 (0.140) 


0.792 (0.775) 



wave functions for both mesons are not singular at short distances and the asymptotic behavior of 
D-state wave function of p is proportional to that of its S-state wave function at small radii. 

With the given wave functions we can obtain the eigenvalues, 6 2 , which in turn give the invariant 
masses of 0.159 GeV for ir and 0.792 GeV for p, while their experimental values are 0.140 GeV and 

0. 775 GeV, respectively. We tabulated the contribution of each term in Eqs.Q-Q in Table [I] The 
spin-spin contribution (< —3^55 >) for pion is larger than the other spin-dependent terms by an order 
of magnitude. However for p the same contribution (< <Pss >) is smaller than the other spin-dependent 
contributions by an order of magnitude. Therefore the spin-dependent term plays crucial role for the 
mass difference between ir and p. The major contribution for pion mass comes from the Darwin (<£_d) 
and spin-spin interaction terms. The magnitudes of these two terms are about same but opposite in 
sign, and therefore they almost cancel each other, resulting in a small pion mass. However for the p, the 
Darwin term is not dominant at all. It is about the same as or less than the off-diagonal contributions 
from tensor (&t) or spin-orbit-tensor (<Psot) terms. Besides that the spin-spin contribution is much 
smaller than the Darwin or the off-diagonal terms, it contributes in same direction as the Darwin term, 
as distinct from the pion case. Overall the p mass is not that small. We can also notice that the kinetic 
energy contributions (< —d 2 /dr 2 + J(J + l)/r 2 >) are almost 1 GeV for ir and 0.5 GeV for p. These 
numbers are much larger than the quark masses of ~ 0.0557 GeV and it justifies to treat the quarks 
relativistically. 

Finally we can conclude that the cancellation between the Darwin and spin-spin interaction terms 
occurs in tt mass but not in p mass and this is the main source of a big difference in two masses. Also 
we can say that our theory is consistent with the chiral symmetry breaking which causes the non-zero 
small mass of pion. 
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